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ABSTRACT 

We study low frequency waves that propagate in a region of layered semi-convection. 
Layered semi-convection is predicted to be present in stellar and planetary interiors 
and can significantly modify the rate of thermal and compositional mixing. We derive 
a series of analytical dispersion relations for plane-parallel layered semi-convection 
in the Boussinesq approximation using a matrix transfer formalism. We find that 
like a continuously stratihed medium, a semi-convective staircase - in which small 
convective regions are separated by sharp density jumps - supports internal gravity 
waves (g-modes). When the wavelength is much longer than the distance between 
semi-convective steps, these behave nearly like g-modes in a continuously stratified 
medium. However, the g-mode period spacing in a semi-convective region is system¬ 
atically smaller than in a continuously stratified medium, and it decreases with de¬ 
creasing mode frequency. When the g-mode wavelength becomes comparable to the 
distance between semi-convective steps, the g-mode frequencies deviate significantly 
from those of a continuously stratified medium (the frequencies are higher). G-modes 
with vertical wavelengths smaller than the distance between semi-convective steps are 
evanescent and do not propagate in the staircase. Thus, there is a lower cutoff fre¬ 
quency for a given horizontal wavenumber. We generalize our results to gravito-inertial 
waves relevant for rapidly rotating stars and planets. Finally, we assess the prospects 
for detecting layered semi-convection using astero/planetary seismology. 


1 INTRODUCTION 


Semi-convection is a fundamental physical process with im¬ 
plications for stellar and planetary evolution and structure. 
In the astrophysical context, semi-convection occurs when 
an unstable thermal (entropy) gradient is stabilized against 
overturning convection by a compositional gradient (e.g. of 


helium) (Schwarzschild & Harm 1958 Kato 1966 Langer 


|et al.|[r983[ |Noels et al.||2010[ ). The requisite physical con¬ 
ditions can occur at the boundary between convective and 
radiative zones in the cores of massive main-sequence stars 
(e.g. B and O stars) ([Langer et al.|[T985| |Langer||1991| and 


helium-burning horizontal branch stars (Castellan! et al. 


1971 Sweigart & Demarque 19721; they can also occur in 


giant planets for which helium and hydrogen are not fully 
miscible ( Stevenson|jl985 Leconte fc Chabrier|2012 Nettel 


mann et al. 20151, and there exists a global compositional 
gradient due to e.g. helium rain out (|Salpeter||1973| [Steven-j 


1975 Wilson & Militzer 20101. For stars, the presence 


of semi-convection has significant ramifications for stellar 
evolution ( |Weiss|1989| jStothers fc Chin|1994| [Gabriel et al.j 
2014); for planets, its presence modifies the global heat con¬ 


tent (entropy) and metal distribution, which might help ex¬ 
plain the inflated radii of hot Jupiters ([Chabrier fc Baraffe 


2007) and the luminosity of Saturn (Leconte & Chabrier 


2013). 


There has been much recent progress in modeling the 


microphysics of semi-convection using both analytical the¬ 


ory and simulations (Rosenblum et al. 2011 Mirouh et al. 


[2012[ [Wood et al.[[2013[ [Zaussinger fc Spruit[[2013[ ). SernI 


convection arises from a double-diffusive instability (Stern 


1960), specifically when heat diffuses much more rapidly 


than composition. Double-diffusive instability can either 
lead to turbulent diffusion or layered semi-convection, de¬ 
pending on the properties of the medium. In this work, 
we focus on layered semi-convection ( [Proctor][1981[ [Spruit[ 
1992 2013 Radko [2003|, which has a staircase structure 


for the density (see Fig. that is a unique signature of 
semi-convection, at least in certain regimes. Layered semi¬ 
convection is analogous to overturning convection but with 
eddies occurring on the scale of the separation between in¬ 
dividual steps in the staircase. The step size separation, 
d, is typically smaller than the scale height, H, with esti¬ 
mates for d/H spa nning orders of magnitude in range from 
10~® < d/H < 1 ( Leccmte fc ChabriCT|[20i2 Zaussinger & 


Spruit [2013 Nettelmann et al. 2015). Semi-convection be¬ 


comes more efficient at tranporting energy with increasing 
step-size, because of the increasing distance convective ed¬ 
dies can travel before re-depositing their energy. 


Layered semi-convection caused by the salt-finger in¬ 
stability has been observed in thermohaline staircases in 
the ocean. However, definitive observational signatures of 
layered semi-convection are lacking in the astrophysical 
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context. Given the high resolution stellar asteroseismology 


data from the Kepler satellite (Gilliland & et. al. 2010 


Christensen-Dalsgaard 20121, as well as Cassini observations 


of density waves excited in Saturn’s rings by internal modes 
of Saturn (Marley fc Porco|1993 Hedman & Nicholson 2013 


Fuller||2014 Fuller et al.||2014 |, it is natural to look for sig¬ 

natures of semi-convection in the oscillation modes of stars 
and planets. 

Our goal in this work is to understand the low frequency 
oscillation modes in layered semi-convection using a local 
plane-parallel approximation. These low frequency modes 
have buoyancy and/or rotation as their primary restoring 
force. We work under the assumption that a semi-convective 
staircase exists, without concerning ourselves with the un¬ 
derlying (and fascinating) physics of how it is generated and 
evolves over time. We find that like a continuous, stably- 
stratified medium, a semi-convective staircase with density 
jumps supports g-modes that are the analogs of internal 
gravity waves. G-modes with wavelengths that are long com¬ 
pared to the step size behave like internal gravity waves, 
though with small corrections that we quantify and which 
may be detectable. However, g-modes with wavelengths of 
order the step size are strongly affected by the discrete na¬ 
ture of the density jumps in the staircase, and their dis¬ 
persion relation differs from that of internal gravity waves. 
Moreover, the discrete nature of the steps in the staircase 
introduces a lower cutoff frequency for g-modes with wave¬ 
lengths smaller than the step size. This cutoff frequency has 
no counterpart in a continously stratified medium. 

We begin by describing the setup of our model for a 
semi-convective staircase in In ^ we perturb around 
this background state and derive analytical expressions for 
the dispersion relation of linear g-modes subject to different 
boundary conditions, under the Boussinesq approximation. 
In we study the effect of rotation on the dispersion rela¬ 
tion of g-modes. In we discuss the potential of detecting 
layered semi-convection using asteroseismology, and we end 
with a summary and a discussion of our results in 


2 PROBLEM SETUP 

We assume a background state that is in hydrostatic equilib¬ 
rium and adopt a plane parallel approximation, so we ignore 
curvature and work on scales that are small compared to the 
radius of the star/planet. Gravity is constant in our setup 
and points in the —z direction so that the equation of hy¬ 
drostatic equilibrium is dP/dz = —gp. 

For the density profile, we assume a semi-convective 
staircase, which has discrete interfaces separating convec¬ 
tive zones. At each interface, the density undergoes a dis¬ 
continuous jump by a value Ap. Since the fluid is convective 
in between any two interfaces (within a stair), the Brunt- 
Vaisala frequency 

__ din Pp-^ 

7P dz dz 

within a stair is vanishingly small compared to the frequen¬ 
cies of interest in the problem, so we set it to zero. 

For simplicity, we shall typically assume that the dis¬ 
tance between adjacent interfaces, d, is constant, though we 
also discuss the effect of relaxing that assumption. We also 



Figure 1. Solid curve — Density profile for layered semi¬ 
convection; dashed curve — Density profile for a stable, contin¬ 
uously stratified medium with Brunt-Vaisala frequency N. Since 
we shall be working on scales that are small compared to the scale 
height, we have used straight line segments to approximate both 
an adiabatic gradient (solid line segments between density jumps) 
and a continuous density profile. 


assume that d H, where H is the characteristic length 
scale for the background equilibrium quantites to change by 
order unity {H is of order the scale height). Under these as¬ 
sumptions, the magnitude of the jump in density, Ap > 0, 
between adjacent stairs is given by 

( 2 ) 

where N is the Brunt-Vaisala frequency given by integrating 
equation ([^ over a single stair, including the interface. Since 
N Ks 0 within a stair due to convection and each interface is 
presumed to be infinitely thin, the only contribution to N is 
from the jump in density across an interface. Also, because 
d <C R by assumption, we generally have Ap/p 1. 

Figure[^shows a schematic of the density profile that we 
consider for layered semi-convection (solid curve) vs. a con¬ 
tinuous density profile with the same value of N (dashed 
line). For the semi-convective profile, the interfaces are 
spaced in intervals of d and the jump at each interface is 
Ap. Between interfaces, p(z) follows an adiabatic gradient, 
due to convection between stairs. 

Having described the equilibrium background state, we 
shall now linearize about it to derive the dispersion relation 
for g-modes. For simplicity, we ignore any corrections to the 
dispersion relation resulting from a finite scale height. Thus, 
we assume the wavelength of the perturbations is much 
smaller than the scale height {X/H <C 1). This limit together 
with the low frequency assumption uj <C Cs/A, where Cs is 
the sound speed, constitute the Boussinesq approximation. 


3 DISPERSION RELATION FOR G-MODES 

We now derive the dispersion relation for the analog of g- 
modes in a semi-convective staircase under the Boussinesq 
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approximation. We initially neglect rotation, but we return 
to the effects of rotation in !|^ Assuming an adiabatic equa¬ 
tion of state and (without loss of generality) motion in the 
X — z plane, the linearized Boussinesq equations are 


0 = 

dSu 
iH ^ 
ddw 

~dr ^ 


dSu dSw 
dx dz 
_ld5P 
p dx 

1 dSP Sp 

- 5 - 9 — 

p oz p 


') dSp dlnPp 
dz 


(3) 


where 5u and 5w are the fluid velocities in the x and 2 
directions, respectively. 

This set of equations can be simplified further by not¬ 
ing that dlxiPp '/dz = 0 within a stair due to convection. 
Thus, the set of equations ([^ within a stair reduces to 


_ dSu dSw 
dx dz 
dSu 1 dSP 

dt p dx 

dSw 1 dSP 

dt p dz ' 


(4) 


from which it immediately follows that \/^SP — 0 

The set of equations @ applies within a stair, but we 
also need to connect adjacent stairs via boundary conditions 
at the interfaces. There are two standard boundary condi¬ 
tions in this case: the contact condition, which states that 
the fluids on either side of the interface must stay in contact, 
and the force balance condition, which states that the forces 
must balance across the interface to avoid infinite accelera¬ 
tions. 

Let us assume for specificity that interfaces are located 
at integer multiples of d, and that the n-th stair lies be¬ 
tween —nd ^ 2 ^ — (n — l)d. Referencing a perturbation in 
the n-th stair by the subscript n, and defining ^ to be the 
displacement vector of a fluid element from its equilibrium 
position, the boundary conditions across an interface are 


^n + l — G (5) 

^Pn + 1 ~ ^Pn T ^nft^P- 


The first boundary condition expresses the requirement that 
the two fluids stay in contact with each other. The second 
boundary condition ensures that the momentum flux is con¬ 
tinuous across the interface and can be obtained by inte¬ 
grating the 2 -component of the momentum equation across 
the interface. All quantities in the equations ([^ are evalu¬ 
ated at the interface located at 2 = —nd, and the boundary 
conditions connect the linear perturbations in the n-th stair 
to those in the (n -|- l)-th stair. 

Next, we would like to express in terms of SPn- To do 
so, we Fourier transform equations in X and t and assume 
normal mode perturbations in the form SP„ oc 
The z-dependence of the eigenfunction follows directly from 
the relation \7^5P = 0. Thus, the form of the pressure per¬ 
turbation for plane waves in the 2 ;-direction can be written 
as 


Combining equation a with the 2 -component of the mo¬ 
mentum equation Q, we see that 

£ri = 

" V " " / 

Here, we have assumed that Ap/p <C 1, in order to be con¬ 
sistent with the Boussinesq approximation, which neglects 
the vertical density variation, except in the buoyancy - gAp. 

The boundary conditions fsf written in terms of A„ and 
Bn are 

An+ie’^^'^ - Bn+ie-'^^'^ ^ An - Bn ( 8 ) 

An + ie^^'^ + Bn + ie-’^^‘^ = An + Bn + - Bn). 

pUJ-^ 

In addition to these equations, which connect the values of 
A and B in adjacent stairs across an interface, we must also 
provide boundary conditions for A and B. Once boundary 
conditions are specified, it is possible to solve for oj and 
derive the dispersion relation. In what follows, we consider 
two different cases: an infinite staircase and a finite staircase 
embedded in a convective zone. 


3.1 Infinite Staircase 


We now derive the dispersion relation for an infinite semi- 
convective staircase. We first consider strictly periodic 
boundary conditions. In other words, we assume that An = 
An.\-m and Bn = Bn+m, so the eigenfunctions repeat every 
m stairs. More general solutions are given in ( 3.1.1| For a pe¬ 
riodicity of m stairs, we see from equations that we will 
in general have a system of 2m equations to solve. Defining 


= gk±Ap/pixd = [N/kSf kxd, 


(9) 


this system of equations can be written in matrix form as 

TP Q 0 ... 0‘ 

0 


0 

Q 0 
P = 


Q 

p 


'vi' 




= 0, v„ = 

Bn 





-1 


-{1 + W^) IT^-l 


, Q = 


( 10 ) 


Here, [ui, ...Vm] is the vector of length 2m composed of the 
coefficients An and Bn and P and Q are 2x2 submatrices. 
P is repeated down the main diagonal of the 2m x 2m ma¬ 
trix m times, and Q is repeated down the diagonal above 
it m — 1 times. There is a corner element, Q, in the lower 
left hand corner, which enforces periodicity, and all other 
matrix elements are zero. 


Equation (101 is not ideal, because the size of the ma¬ 


trix grows without bound as m increases. However, one can 
write down the solution of block tridiagonal systems in the 


form (101 by making use of a recurrence relation ( Molinari 
|2008[ ). For instance, we can write {An+i, Bn+i} in terms of 
{AnjBn} as 


( 11 ) 




An + l 

Bn 


Bn+\ 
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4 Mikhail A. Belyaev, Eliot Quataert and Jim Fuller 


where the transfer matrix T is defined by 

'^1 _l_ e"*" iii/2,-fc|d 

f 1 W‘‘ \ „k I d 


T = 




( 12 ) 


The system of equations ( |10[ | with periodic boundary con¬ 
ditions can be written succintly using the transfer matrix 
as 


(13) 




An 

Bn 


Bn 


In order for equation (131 to have a non-trivial solution, we 
must have 


det(T’" - /) = 0, 

where I is the 2x2 identity matrix. 


(14) 


One way to solve equation (141 is to compute — I 
explicitly and take the determinant. Computing T™ can be 
accomplished by diagonalizing T, i.e. by writing the transfer 
matrix in the form T = UDU~^, where D is a diagonal 
matrix. The transfer matrix raised to the power of m is then 
simply given by T™ = UD’^U~^. Alternatively, equation 
(141 can be solved by noting that it implies both of the 
eigenvalues of the matrix T™ are equal to 1, i.e. A[T™'] = 1 
for both eigenvalues. This is equivalent to the expression 
(A[T])’" = 1 for the eigenvalues of the transfer matrix T. 

Solving equation (|14[), we find after some algebra 


/ W 

1 = ( cosh(fcj^d)-^ sinh(fc^d) ± 


sinh(fc^d) 


11/4 


-b 1 J sinh(fe^ci) — cosh{k±d) 


(15) 


We can g et rid of the functional dependence on k±d in 
equation (151 by assuming solutions of the form = 
2coth(fc^d) — 2Ccsch{k±d), where C is a constant to be 


solved for. Making this substitution, equation (151 simpli¬ 
fies to 


(C±y/C^ - l)"* = 1- 


(16) 


Equation (161 will be satisfied if C ± y/C^ — 1 = 


where n is an integer. Solving this equation for C we have 
C = cos(27rn/m). (17) 

Note that since cos(a;) = cos(—i) there are only m/2-|-l 
distinct values of C, where we round down if m is odd. We 
should in fact choose 0 ^ n ^ m/2, because if n = m/2 the 
eigenfunction repeats every other stair, and higher values of 
m are “aliased” to lower values. This is analogous to what 
happens in signal processing when frequencies higher than 
half the sampling frequency are aliased to lower frequencies 
(i.e. the Nyquist-Shannon sampling theorem). 

Given equation ( |17[ ), the dispersion relation for periodic 
solutions is 

k I d 


(18) 

^2coth(fc^d) — 2 cos(27rn/m) csch(fcxd) J ’ 

where iV^, defined in equation ||^ , is the background Brunt- 
Vaisala frequency that would correspond to a smooth den¬ 
sity gradient in the absence of semi-convection. 


3.1.1 General Dispersion Relation for the Infinite 
Staircase 


The solution in equation (181 came from requiring that 


the eigenfunctions repeat every m stairs. This strict peri¬ 
odicity requirement led to a discrete spectrum of modes 
with a finite number of frequencies. We now relax the as¬ 
sumption of strict periodicity and look for solutions with 
[An+m,Bn+m] = [An, Bn]. The dispersion relation in 
this case is given by the expression 

det(T’" - = 0, 


(19) 


where in general if is a complex constant, ip = ipr + iipi with 
ipr restricted to the range Q ^ ipr < 1 (without loss of gen¬ 
erality). Physically, ipr represents the phase shift a solution 
undergoes over m stairs. Strict periodicity requires ip = Q 


(equation (14|), but we can still have purely oscillatory so¬ 


lutions under the less stringent requirement oi ipi = 0. So¬ 
lutions with ipi ^ 0 are spatially growing/decaying and are 
thus analogous to evanescent waves. 

Following the derivation given in §3.11 and again as¬ 
suming solutions for W in the form W — 2coth(k±d) — 
2C csch(fcxd), where C is a constant, the dispersion relation 
(191 reduces to 


C ± yJC^ - 1 

g27ri'i/’ / m 


= 1 , 


( 20 ) 


which is the analog of equation ( |16[ ). Taking C = cos(6'), 
where 9 is in general complex, equation (201 reduces to 


= 1 , 


which has solutions 


2'K{n±ipr) . 2'Kiipi 
a = -±-, 


( 21 ) 


( 22 ) 


where n is an integer. 

Considering just oscillatory solutions for the time being, 
which have ipi = 0, we see that the spectrum of modes for the 
infinite staircase is actually continuous, and the discreteness 
obtained in the previous section is an artifact of imposing a 
periodicity constraint. Thus, the dispersion relation for the 
infinite staircase can be written as 

k±d 


2 coth(fcxd) — 2 cos{kzd) csch(A:xd) 


(23) 


where we have substituted 9 — kzd. The reason for this 
substitution will become clear in the next section when we 
treat kz as the vertical wavevector for a g-mode propagating 
through the staircase. For purely oscillatory solutions, we 
have the requirement that k^ is real. 

3.1.2 Analysis of Infinite Staircase Solution 

In order to better understand the frequency shifts induced 
by the presence of a semi-convective staircase compared to a 
continuous stratificati on, it is instructive to take the contin¬ 
uous limit of equation (231. Taking the limit y/k^ + kid —0 
we can write 




k/ 


kl+kl 


+ O [{kl + kl)d^ 


(24) 


The leading order term, which we shall refer to as cjq, gives 
the dispersion relation for Boussinesq gravity waves in a 
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Figure 2. Left: Exact and approximate (series expansion from equation l |25[ l) curves for Alj/N as a function of kzd for k±d = {0 .5,1,1.5}. 
Right: Plot of of lu/N as a function of kz/k± for various values of k^d. The continuous limit {k^d = 0) given by equation | |24^ is shown 
by the blue line. 


continuously stratified medium. Thus, the dispersion rela¬ 
tion for the semi-convective staircase converges to that of a 
continuous stratification when the wavelength of a mode is 
much larger than the distance between stairs. 

We can ask how much does the frequency of an eigen¬ 
mode for the semi-convective staircase deviate relative to 
that of a Boussinesq gravity wave with the same values of 
fcy and fcz? The most straightforward quantity to consider 
is Aoj = uj — loq, which is the deviation in frequency of the 
staircase from that of a continuously stratified medium. Ex¬ 
panding uj in powers of k^d and k^d, we can approximate 
the fractional change in frequency as 


tion from the continuous case increases with k^d, for a given 




(25) 


This expression is good to a few percent for values of k±d < 1 
and for any value of kzd. The latter is true because kzd ^ n, 


due to aliasing () 3.11, and higher powers of kzd are damped 
by numerical prefactors for kzd ^ tt. For kzd < 1, the fourth 
order term in kzd can be neglected, but it is needed to pro¬ 
vide a good approximation in the range f ^ kzd < tt. 

The left panel of Fig compares the exact value of 
Auj/N as a function of kzd for k±d — {0.5,1,1.5} (solid 
lines) with the series approximation from equation (251 
(dashed lines). Even for k±d = 1, equation (251 is accurate 
to within a few percent for Auj/N across the entire range 
of values for kz. Notice that Alj > 0 always, which can also 
be seen from equation ( |25| l and implies that the staircase is 
“stiffer” than a continuously stratified medium. 

The right panel of Fig compares uj/N as a function 
of kz/k±^ for various values of k±d. For a continuous density 
stratification, ljj/N = ujq/N is a single curve given by equa¬ 
tion (241. For a stratified staircase, however, the frequencies 


are higher than for a continuous density stratification. The 
black curves corresponding to the staircase dispersion rela¬ 
tion in the right panel of Fig [^terminate abruptly, because 
the maximum value of fcj/fcy is vr/feyd as a result of the 
condition kzd tt. 

From the right panel of Fig it is clear that the devia¬ 


k±d, which is predicted by equation (251. Also, there is a 


minimum value of uj, which occurs at kzd = tt. For a given 
fcyd, this cutoff frequency is given by 


cUc k±dsinh{k±d) 
2(cosh(fcyd) + 1) 


(26) 


In the limit k±d 1, the cutoff frequency is ujc ~ Nk±d/2. 
In a continuously stratified medium, this would correspond 
to the frequency of an internal gravity wave which has a 
vertical wavelength of order the separation between stairs 
and a horizontal wavevector equal to k±. 

Although oscillatory solutions below the cutoff fre¬ 
quency, uJc, are not allowed, it is possible to have evanescent 
solutions. These correspond to a complex vertical wavevec¬ 
tor, kz. For instance, if we take kzd = tt -\- iad, where 
—oo < a < oo, equation (|23|) becomes 


2 _ 7vr2 

u = N 


k±d 


2 coth(fcyd) -|- 2 cosh(ad) csch{k±d) 


(27) 


Thus, we can have any value of the frequency in the range 
0 < u) < cJc, by choosing the magnitude of a. The sign 
of a doesn’t matter for the dispersion relation, but it does 
affect whether the waves grow or damp with increasing z. 


From equations (191 and (221, we see that the exponential 


growth/damping rate per unit length for an evanescent wave 
is simply 1/ad. The interpretation of these evanescent so¬ 
lutions with complex kz is that if a wave with frequency 
cu < ojc is incident on a semi-convective region it will be 
reflected. 

Another point worth mentioning about the right panel 
of Figj^is that the curves for ur/N flatten out for k±d > 1 
and become nearly constant as a function of kz. Indeed, in 
the limit of 1, the exact dispersion relation ( |18[ ) is 

independent of kz 


2 N^k±d gk±Ap 


(28) 


2 2p 

This is the standard dispersion relation for an interfacial 
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gravity wave in the Boussinesq approximation {/S.plp <C 1) 
in an infinitely extended fiuid. This makes sense, because the 
limit k±d^ 1 corresponds to the case when the wavelength 
is much shorter than the distance between stairs, so the wave 
dynamics is set by a single interface. 


3.2 Finite Staircase Embedded in a Convective 
Medium 

We now derive the dispersion relation for a semi-convective 
staircase of finite vertical extent with m interfaces, embed¬ 
ded in a convective medium. We again need to solve a system 


of equations similar to ( 101 , but with boundary conditions 


Ao = 0 at the top of the staircase and Bm = 0 at the bottom. 
These boundary conditions come from requiring the pertur¬ 
bations to decay as z ^ ±oo, since the convective medium in 
which the staircase is embedded does not support g-modes. 

Using the transfer matrix T defined in equation (111 
and applying boundary conditions, we can write 


(29) 


■ 0 ■ 





0 


In order for equation (291 to be true in general, the lower 


right hand corner of the 2 x 2 matrix T™ must equal zero. 
Thus, the dispersion relation becomes 


= 0 , 


(30) 


where our row and column indices take the values of 0 and 1 . 
Computing T’" explicitly using diagonalization and making 
the substitution 


— 2 coth(fcj^d) — 2cos{9) csch{k±d), 
the dispersion relation can be written as 

e*"*® {cos(6) coth(fcj_d) — csch(fe_Lci) -|- isin(0)) = 
^-imB coth(k±d) — csch(k±d) — . 


= 0, 


Equation ( |32[ | is of the form 

^imO ^ — imO 

so we can simplify it to read 

sin(m 0 -(-<())= 0 , 


(31) 


(32) 


(33) 


(34) 


where 0 ^ (j> < n is implicitly defined through equations 


(321 and (331. We can convert equation (341 into an explicit 


expression for 9 by using the trigonometric angle addition 
formula. Assuming (f) J 0 and 0 7 ^ 0, we can write 

cot(m^) -I- cot 0 coth(A:j_d) — csc9 csc]i{k±d) = 0, (35) 

where we have used the relation cot{(j>) = cot 9 coih{k±d) — 
CSC 9 csch(fcj_d). 

We can turn equation ( |35[ ) into a polynomial expression 
for cos{9) by using the relations 


cos{m9) = Tm{cos9) 
sm{m9) = Um-i{cos9) sinS, 


(36) 


where Tm and Um are Chebyshev polynomials of the first 
and second kind, respectively. Using these identities, equa¬ 
tion (|35[l can be written as a polynomial in cos{9) 


Tm(cos9) -I- 

(cos 6 coth(fcj^d) — csch(fc^d))[/m-i(cos 6 ') = 0. (37) 


The dispersion relation is given by equation ( |31[ ) , where 
cos(9) takes the values of the roots of the polynomial equa¬ 
tion (37l. Since equation (371 is an m-th degree polynomial, 
there will in general be m branches of the dispersion relation 
for an embedded staircase with m interfaces. 

An alternative to equation (371 is to use equation (341 
to solve for 9 directly. In particular, we can write 


Txn (j) 

m m 


(38) 


Taking the limits of m ^ 1 (large number of steps) and n ^ 
1 (small vertical wavelength relative to size of staircase), we 
see that 9 = 'Knim -|- since 0 ^ ((> < tt. In this 

limit, cos 9 ~ cos( 7 rn/m), and the dispersion relation for a 
finite staircase with m interfaces converges to the dispersion 
relation for an infinitely extended staircase (equation ( |18[ )), 
assuming a periodicity of 2 m stairs. 

The convergence to the infinite staircase is not surpris¬ 
ing, because in the limit 1 many vertical wavelengths fit 
between the upper and lower boundaries, and their frequen¬ 
cies are only weakly influenced by the boundary conditions. 
However, one might wonder why cos 9 —>■ cos(7rn/m) rather 
than cos(27rn/m) (i.e. why the finite staircase with m 3> 1 
and n 1 is equivalent to a staircase of periodicity 2m 
rather than m). The reason is that the finite staircase allows 
“antiperiodic” solutions as well as strictly periodic ones. An 
antiperiodic solution of m stairs obeys the equation 


An 


An 

Bn 

1 

1 

Bn 


(39) 


The analog of equation (131 that permits both periodic and 
antiperiodic solutions is 


det(T’" ± /) = 0. (40) 


The solution to this equation is 

= 2 coth(fc^d) — 2 cos( 7 m/m) csch(A:^d), (41) 

which contains a factor of cos(7rn/m) rather than 
cos(27rn/m). 


3.3 Nonuniform Step Sizes 


Up to now we have assumed that the step-size, d, is a con¬ 
stant. We now relax that assumption and ask what happens 
when the step size is nonuniform? Assuming periodic solu¬ 
tions, we find that uj is an even function of k±d, where d 
is a characteristic value of the step size derived by an ap¬ 
propriate average over the the step size distribution. Thus, 
for k±d —7 0 , the frequency deviation from the continuous 
case must be at least second order in k±d, which we showed 
explicitly for the case of equal step sizes (equation (HI})- 
To treat the problem of unequal step-sizes, we can still 
use the matrix transfer formalism. In particular, if the step 


size of the n-th stair is equation (14 1 can be written as 


det(TiT 2 ...Tm — J) — 0, 


(42) 


where the transfer matrix T„ is the same as the transfer 


matrix T from equation (111, but with d replaced by d„. It 


is also useful to define = a„d where the a„ are numerical 
coefficients. Thus, the n-th transfer matrix can be written 
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= 




' 2 ^ 






(43) 


where 6„ = exp(a„fc^d) and = an{N/io)‘^k±d. 

Equation (421 is an implicit equation for 
uj{k±d,ai, ...,am)- To show that uj is an even function 
of k±d, with the a„ coefficients fixed, we need to demon¬ 
strate that oj{k±d) = ij{—k±d). Since uj{k±d) satisfies 
equa tion (421, it is clear that uj{—k_i_d) also satisfies equa¬ 
tion (421, but with hr^-^ l/&n and 


—W^. Thus, the 


analog of equation (421 for uj{—k±d,ai,...,am) is 


det{TlT*2...Tl - I) = 0, 


where 


rri* _ 

J. yy 


i^-9) 


' ‘2bn 






(44) 


(45) 


We see that TJ) is the same as Tn but with the indices 
on both the row and columns flipped (i.e. T^c = Tt-c,i-r)- 
From the definition of the 2x2 determinant, it is clear 


that equations (421 and (441 are equ ivalent. Consequently, 
uj{—k±d) also satisfies equation (42 1 , which concludes the 


proof that uj must be an even function of k±d. 

For a staircase with a hnite number of steps, one may 
also wonder whether going to a solution with nonuniform 
step size, but the same number of steps, changes the num¬ 
ber of solutions for cu that satisfy the dispersion relation. Re¬ 
gardless of the boundary conditions, the maximum number 
of distinct solutions for u!^ (or alternatively for cu with the 
assumption cu ^ 0) is equal to m, the number of steps. This 
is because multiplication of an input vector by a set of m 
transfer matrices effectively generates an m-th degree poly¬ 
nomial for cu^, since each transfer matrix is linear in 1/cu^. 
The dispersion relation for a finite staircase embedded in a 


convective medium has m distinct roots (e.g. equation (34|), 


which means it is complete. Going to non-uniform step sizes 
would shift the positions of these roots, but the number of 
solutions would remain the same. This conclusion explicitly 
demonstrates that non-uniform step sizes cannot produce 
“splitting” of mode frequencies. 


Fourier transforming in time (i.e. d/dt —>■ —icu) and 
changing variables to the fluid displacement vector, de¬ 
fined via 


X c 

= - = -zcuC, 


equations (461 can be written as 
V ■ ^ = 0 


-cu^^ = -b2icu(r2 X ^). 


Taking the curl of equation (491 we have 


(47) 

(48) 

(49) 

(50) 

(51) 

Substituting equation ( |50| l into ( |51| l, we hnally arrive at 

cuV^^ = 4(n ■ V)(n ■ V)^ (52) 

Assuming solutions in the form 

^ = f{z)exp[i{kj_ ■ X - ujt)], (53) 

where fen is the wavevector in the x — y plane (i.e. the plane 
perpendicular to the stratification) and substituting equa¬ 
tion (531 into equation (511, we have the following second 


—w^(V X ^) = ■ V)^, 

and taking the curl again we have 

= -2iuj{n ■ V)(V X ^). 


order ODE for f{z 

dz^ 


{ -k\ + ) f{z) = 


4 ( ik±Q,± + Dj 


dz 


ik±fl±+Q.^-^] f{z), (54) 


where Dx = ft ■ k±/k±. 

The solutions to this ODE have the form f{z) oc 
exp{kzz), where k^ is defined as 


kz = k± 


4ifl±flz ± — 4 (S2^ -|- fl|) 




(55) 


In what follows, we restrict our attention to ft parallel to z 
and ft perpendicular to z and focus on the infinite staircase. 


4 ROTATION EFFECTS 


We now address how to include the effect of a constant ro¬ 
tation rate in the dispersion relation. We assume that the 
rotation rate is far below breakup and so the centrifugal 
force can be ignored. In this limit, the Coriolis force is the 
key new restoring force and introduces inertial waves into 
the problem. The linearized fluid equations inside a stair in 
the corotating frame are 


V • = 0 


dSv 

ILL 


VSP 

P 


— 2ft X 5v. 


(46) 


The two major differences with respect to equations (|^ are 
the inclusion of the Coriolis force and that we can no longer 
assume the motion is in a 2D plane. The latter is true if ft 
is misaligned with the direction of density stratification. 


4.1 Rotation Axis Parallel to « 


As we show in Append ix |A1[ the boundary conditions be¬ 
tween stairs, equations (|5|and 0 . are unchanged for a ro¬ 
tation axis aligned with « if we replace fex with fc: 


k = k± 



(56) 


Thus, for an infinite staircase, the dispersion relation can be 
written as 


2 10-2 
oj = N 


kd 


2 coth(fcd) — 2 cos{kzd) csch(fcd) 


(57) 


Because equation J 53 has the same form as equation 
|, the results of ^3.1.2| carry over. For instance, the analog 
of equation (241 in the continuous limit, kd <C 1, is 


-I- 


+ O [{kdf 


(58) 
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k±d = 0.5, kjk±=0 k±d = 0.5, kjk±=.5 




Figure 3. Comparison of exact and approximate solutions for ui/N using the exact dispersion relation j57[ l and an inverse RMS weighting 
of approximations l|64[l and l|65||. The left panel shows fcz/fcx = 0, and the right panel shows kz/k± = 0.5. 


Substituting for k and using the definitions of k±, k^, and 
cjo from §3.1.2[ we have 


2 AT^ 

luq = -ly 


wg + k^(LU^ - 4^1) ■ 


Solving for tuo yields 


2 _ fcliv" + 4k^ni 

UIq - 


kj + 


(59) 


(60) 


The exact Boussinesq dispersion relation for plane waves in a 
continuously stratified medium and an arbitrary orientation 
of fl with respect to the direction of stratification is 


2 fci7V" + 4(fcj 

LJ — - 




+fc 2 


(61) 


Comparing this with equation (601, we see that the disper¬ 
sion relation (571 correctly recovers the continuous limit. 

In addition to computing wo, we can also compute the 
leading order approximation for the effective wavenumber 


(62) 


by plugging equation (601 into equation ( |56[ |: 

l + {2^lz/Nf{kz/kj_f- 


1 - {2nz/NJ 


We see that k^ ^ with equality only if = 0 and that 
fco diverges for iV^ = 

Next, we compute the leading order correction to the 
dispersion relation in the limit kd —>■ 0. Expanding cj to 
second order in kd for fed <C 1, we have 


uj _ I fc2 1 {kodf 
W “ V k'^ + kl ^ 24 ojo/N ’ 


(63) 


where cuo is given by equation ( |60[ ). Note that we cannot 
replace the first term on the right hand side of equation (631 
with cJo, because k contains deviations from feo (equation 
at second order in k±d and ultimately we want an 
expression in terms of k±d not kd. Solving for Atj/tuo gives 


Alu ^ I / {k±df \ 
EJ 24 \1- {2nz/N)^) 


(64) 


Note that for N too is finite but Ao; —>■ 0. This limit 

represents purely inertial waves unaffected by the staircase. 


The ma jor modification to Au in equation (641 relative 
to equation (251 is the 1 — {2Q,z/NJ term in the denomina¬ 
tor, which makes Acu/tuo blow up for N’^ = However, 
ko also blows up for = 40,^ (equation (621), so our ex¬ 
pansion in fed ^ 1 is not valid. To understand what happens 
when iV^ = we can set u) = N + Ato, since u>o = N 


is an exact solution to equation (601 in this case. Thus, Aco 


represents the difference in frequency between the staircase 


and a continuously stratified medium. From equation (571, 
one finds that for fed ^ 1 


Alj 

w 


kid 


4^3 \ ^l + {kz/kM^ 


nEo = 40.1, (65) 


and in particular, there is no divergence in Atj. Addition¬ 
ally, there is no divergence in fc for iV^ = 4fl^, and fe can 
be estimated by plugging u) = N + Auj with At*; given by 
equatio n (|65[ l into equation ( |56[ ). 

Fig 1^ shows the exact value of uj and the approximation 
LO = wq + Ao; for the staircase as a function of Q^/N. For this 
figure, Alj is computed using an inverse root mean square 


weighting that combines approximations (641 and (651. 


4.2 Rotation Axis Perpendicular to 2 

In the case when the spin axis is perpendicular to 2 (i.e. in 
the X — y plane), we can assume without loss of generality 
that O = O^x. The jump conditions at an interface are 
now different from those given in equations § and (|^, and 
we derive the jump conditions for ft perpendicular to z in 
Appendix I A2| Using these jump conditions and defining 

/ \Jlj‘^ - 401 {kx/k±f 
“ ^ 1 UJ 
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the dispersion relation for Q, = Q.j:X can be written as 

^ ^ kd 

/ V 2 coth(feci) — 2 cos(fezd) csch(fcd) 

Equation (1671) can be rearranged to read 


uj^ = nU% 


(67) 


uj^ -4klli ^] = 




kd 


2 coth(fcd) — 2 cos(fcz(i) csch{kd) 


( 68 ) 


We can check that in the continuous limit, expression 
( |68[ ) reproduces the dispersion relation for a continuous den¬ 
sity stratification with a rotation axis perpendicular to a. 
Taking fed —>■ 0, we have 


2 _ ,o2 / kx 

u — 4Qx I ^2 


-b W 


k^ -b 


+ o [{kdf 


(69) 


Expanding k according to equation (|66 
duces to 


equation (69) re- 


2 

UJq 


N^k 


2'^2 


G. -b 


(70) 


Upon comparison with equation (611, it is clear that the dis¬ 


persion relation (671 correctly recovers the continuous limit. 


We can also solve for k to leading order 

{2nx/Nf{kx/k^f (1 -b {kGk±f) 


f 2 7 2 

ko — 


1 - 


1-b (2n„/iV)2(fc,/fcx)2 


(71) 


from which we see that fcg SS k\ with strict equality only if 

klxkx — 0 . 

Next, we solve for the leading order corrections to u) 
and find that for fed —> 0 


Ao; 


(fci + kl)d^ 


tuo 24 \^[i + (2n^/jv)2(fc,„/fc^)2] 


(72) 


In contrast to the case when the spin axis is aligned with z, 
there is no divergence of Atu/tJo, since fco does not diverge 


(equation (711). Also notice from equations (72) and (25) 


that Atu/wo is always smaller (for the same values of k^d 
and kzd) for a spin axis perpendicular to 2 compared to the 
case without rotation. Fig shows Alu/ui as a function of 
klxkx/Nk± for kz/k± = 1 and k±d = .5. The approximate 
and exact solutions are nearly on top of each other. 


5 APPLICATION TO ASTEROSEISMOLOGY 


The previous sections show that the presence of lay¬ 
ered semi-convection introduces quantifiable wavelength- 
dependent changes to the dispersion relation for g-modes 
compared to a continuously stratihed medium. However, in 
order to find an unambiguous signature of semi-convection 
in asteroseismic (or planetary seismic) data, it is neces¬ 
sary to incorporate our results into a realistic calculation 
of a stellar/planetary mode oscillation spectrum. A large 
a priori uncertainty in such an analysis is the number of 
density jumps per scale height in layered semi- convection, 
H/d, which can be anywhere from a few to 10® (Leconte & 


Chabrier|2012||Zaussinger & Spruit|20131 |Nettelmann et al. 

20151). 


k±d = Q.5, kjk^ = \ 



Figure 4. Fractional frequency difference, Auj/ui, between a stair¬ 
case with rotation perpendicular to stratification and a continu¬ 
ously stratified medium for k±d = .5 and kz/kj_ = 1. There is 
excellent agreement between the exact dispersion relation, equa¬ 
tion l|67[), and the approximate solution, equation \72\. 


When H/d is of order unity, one can incorporate the ma¬ 
trix transfer formalism into global seismology calculations 
directly by computing the jump in the eigenfunction across 
each step or by modeling the steps as sharp but continuous 
increases in the density and directly integrating the differ¬ 
ential equations across them. However, modeling the jumps 
individually is likely to become prohibitively expensive as 
H/d becomes large. In that case, ^ and show that the 
layered semi-convective medium has a modified dispersion 
relation compared to a continuously stratified medium. 

A unique signature of g-modes in asteroseismology is 
typically that adjacent global normal modes have a nearly 
fixed period spacing. However, there are corrections due to 
e.g. a finite scale height and a cavity size that depends on ra¬ 
dial wavenumber. Here, we assess the effect of layered semi¬ 
convection on the period spacing relation of g-modes com¬ 
pared to a continuously stratified medium, in order to de¬ 
duce the seismological signatures of layered semi-convection. 

In what follows, we consider a cavity of fixed size with 
reflecting walls in a plane parallel medium that has constant 
gravity and Brunt-Vaisala frequency. We then ask how the 
constant period spacing relation is affected by layered semi¬ 
convection within the cavity? For our first example, we con¬ 
sider a g-mode cavity that has layered semi-convection ev¬ 
erywhere (left panel of Fig. [^. This is what one might expect 
for a giant planet with semi-convection throughout the in¬ 
terior (or at least throughout the stably stratified region). 
For our second example, we consider a cavity that is partly 
stably-stratified and partly semi-convective (right panel of 
Fig.||). This is what one would expect for e.g. a star under¬ 
going helium burning with a convective core. For simplicity, 
we do not consider rotation, but it is possible to include it 
in the analysis. 

For a fully semi-convective cavity with constant N, kz 
is constant for a given value of the period P = 27r/cu. Be¬ 
cause the cavity has reflecting walls, an integer number of 
wavelengths must fit inside the cavity. Consequently, 

kz = 2'Kn/L, (73) 
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fully semi-convective partly semi-convective 



0 L 0 L, L 


z z 

Figure 5. Schematic of a g-mode of a given frequency trapped between bounding walls (gray) in a cavity of size L (in the z-dimension). 
Left panel: the g-mode cavity consists of a single zone of layered semi-convection with constant N. Right panel: the g-mode cavity 
contains a zone of layered semi-convection of size La adjacent to a stably stratified zone of size L — La- The value of the Brunt-Vaisala 
frequency in the stably stratified zone is also N. However, due to a modified (stiffer) dispersion relation in the semi-convective zone, 
is greater there than in the stably stratified zone for a given ui. 


where n (S 1 is an integer and L is the size of the cavity. As 
we have shown in (3.1.2 in the limit d —>■ 0, we recover (to 


zeroth order in d) the dispersion relation for a continnously 
stratified medium with Brunt-Vaisala frequency N. Written 
in terms of P and n, the dispersion relation for a continuosly 
stably stratified medium is 


(27r)^n 

Nk±L' 


(74) 


The period spacing relation can be obtained by differentiat¬ 


ing equation (741 with respect to n, which gives 


APo = 


Nk,L' 


(75) 


This is the period change in going from n to n -|- 1 nodes in 
the z dimension. Note that APq is a constant independent 
of n (i.e. the mode frequency). 


From equation (251, it is clear that for a semi-convective 


staircase there is a second order correction to the frequency, 
and hence to the period, in both fcyd and k^d compared 
to a continuously stably-stratified medium. Astrophysically 
relevant g-modes typically have spherical harmonic degree 
I — 1 OT I = 2. Since k± ~ i/r* in our setup, k±d ~ Idjr^. 
Because a semi-convective staircase extends over a region of 
the star that is a small fraction of the stellar radius, the pe¬ 
riod spacing correction in the semi-convective staircase due 
to feyd is too small to be observable. However, k^ is larger 
than k±_ by a factor of {n/l){r,/L), and for n ^ 1 the pe¬ 
riod spacing correction due to semi-convection is potentially 
observable. We work out that correction below: first for a g- 
mode cavity that has layered semi-convection throughout 
(left panel of Fig. and subsequently for a cavity that is 
partly semi-convective (right panel of Fig. [^. 


For layered semi-convection, we can use equation (251 
to write 


2'Kkz 

iVfc, 




(76) 


which is valid to second order in d in the limit kzd <C 1. 
Equation (731 still applies, because by assumption the prop¬ 
erties of the medium are independent of 2 within the cavity, 
so we can write 


(271 


Nk±L 


1 / 2'Knd\J 
^~2A\~ ) 


(77) 


Differentiating equation (771 with respect to n, the pe¬ 
riod spacing relation with the correction for layered semi¬ 
convection is 


AP 




(78) 


From equation ( |78[ ), we see that the g-mode period 
spacing is smaller for the staircase than for a continuously 
stratified medium. We also see that the g-mode period spac¬ 
ing decreases with increasing n, or alternatively decreasing 
frequency. Moreover, the correction due to layered semi¬ 
convection becomes more pronounced as n increases, which 
is the opposite of e.g. the correction due to a finite scale 
height, which becomes smaller with increasing n. Although 
our approximation ( |78| l is not valid for 2'Knd/L ~ 1 (or 
equivalently k^d ~ 1), we can understand what happens in 
this limit by looking at the right panel of Fig[^ From the 
figure, we see that at the cutoff frequency, w = uic, the quan¬ 
tity duj/dkz = 0. This implies that dP/dn = 0 at a; = ojc, 
so that the period spacing tends to zero as we approach the 


© 0000 RAS, MNRAS 000, 000-000 























































G-modes in Layered Semi-Convection 11 


cutoff frequency. Below the cutoff frequency waves do not 
propagate in the layered semi-convective medium. 

Next we consider the period spacing relation for a g- 
mode cavity which contains a region that is stably strati¬ 
fied adjacent to a region undergoing layered semi-convection 
(right panel of Fig. [^. We reiterate that this is common in 
post main sequence stellar interiors. One can show that for 


kzd 1, the analog of equation (781 for a g-mode cavity 
that is partly semi-convective is 


AP = APo 1 


1 / Tind'C Ls 

~ 2 ) T 


(79) 


where Ls/L is the ratio of the size of the semi-convective 
region to the total size of the g-mode cavity. Once again, we 
we see that the period spacing is smaller for the model with 
layered semi-convection and AP decreases for larger n, i.e. 
lower frequency g-modes. 


Comparing equations (781 and (791, one sees that the 


magnitude of the correction to the period spacing relation 
is proportional to Ls/L, so the correction becomes larger 
as semi-convection fills a fractionally larger volume of the 
cavity. Another major difference compared to a fully semi- 
convective cavity is that modes with frequencies lo < ujc can 
propagate within the cavity, but only in the stably strat¬ 
ified region, not in the region of layered semi-convection 


(see equation (26l). Thus, there is an approximately con¬ 


stant period spacing AP ~ {2n)'^/Nk±L for oj ^ ujc and 
a different approximately constant period spacing AP ~ 
/Nk±{L — Ls) for a; ^ uic- The transition between the 
two period spacings occurs at frequencies w ^ Wc (equation 

(HH). 


6 DISCUSSION 

We have analyzed the dispersion relation of g-modes (buoy¬ 
ancy waves) propagating through a plane-parallel fluid 
which has the structure of a layered semi-convective stair¬ 
case. We assumed that the size of the stairs in the semi- 
convective staircase was small compared to the scale height. 
The fluid was assumed to be perfectly convective inside of 
a stair (vanishing entropy gradient and Brunt-Vaisala fre¬ 
quency), and we imposed a small (relative to the background 
density), discontinuous jump in the density at the interfaces 
between stairs (see Fig. [^. Although this simplified model 
abstracted away the rich, double-diffusive physics leading to 
semi-convection, its simplicity allowed us to derive analytical 
dispersion relations for g-modes in a semi-convective stair¬ 
case. Such analytical results are particularly relevant given 
the rich asteroseismic data emerging from e.g. the Kepler 
mission and can serve as guide posts for understanding the 
effect of a region of layered semi-convection on the propaga¬ 
tion of g-modes in a star or planet. 

Using a matrix transfer formalism ( Molina^|2008 (, we 
were able to find analytical solutions for the dispersion rela¬ 
tion of g-modes in layered semi-convection under either the 
assumption of an infinitely extended staircase (( 3.1|3.1~ I 
or a staircase of finite vertical extent embedded in a fully- 


convective medium ()3.2|. We were also able to quantify 


the effect of rotation on the dispersion relation of g-modes 
and explicitly determine how the properties transition from 
pure buoyancy waves to inertial waves as the rotation rate 


increases (@- Note that the latter are unaffected by the 
presence of a staircase. 

The dispersion relation for g-modes in layered semi¬ 
convection is stiffer - i.e. the frequency is higher for a given 
wavelength - than for internal gravity waves in a continu¬ 
ously stratified medium. Modes with wavelengths that are 
long compared to the step size behave like gravity modes 
with a leading order fractional correction to the frequency 
of 0[{d/Xf'], where A is the wavelength and d is the dis¬ 
tance between layers in the staircase (see equation (25l). 


On the other hand, frequencies of g-modes with wavelengths 
comparable to the step size are strongly affected by the dis¬ 
creteness of the density jumps in the staircase. One of the 
most dramatic manifestations of this difference is a lower 
cutoff frequency for the staircase, which does not exist for a 
continuously stratified medium. The effect of this cutoff fre¬ 
quency is to prevent the propagation of modes with vertical 
wavelengths shorter than the step size in the staircase. 

The stiffer dispersion relation and the presence of a cut¬ 
off frequency both have implications for detecting layered 
semi-convection using asteroseismology (see In particu¬ 
lar, layered semi-convection decreases the period spacing of 
g-modes relative to that expected for a continuously strat¬ 
ified medium (il and equations |78| and |79[ ). Moreover, the 
decrease in the period spacing is frequency dependent, be¬ 
coming larger for lower frequency (shorter wavelength) g- 
modes that are more affected by the discrete density jumps 
in the semi-convective staircase. 

The change in g-mode period spacing due to layered 
semi-convection we predict in equations ( 78|79 1 can be dif¬ 
ferentiated from other physical effects. Non-WKB correc¬ 
tions to the period spacing of g-modes become smaller for 
lower frequency waves, in contrast to the semi-convective 
corrections which become larger for lower frequency waves. 
In addition, mode trapping (which can also be created by 
composition gradients but is due to non-WKB effects on 
the mode frequencies) creates localized “dips” in the period 
spacing (see, e.g., |Papics et ah (2015l). Rapid rotation often 
creates a period spacing that changes nearly linearly with 
radial order n (Bouabid et al. poT^ , with an opposite slope 
for prograde and retrograde modes. 

Several well-measured period spacings in pulsating stars 
with convective cores (which may contain semi-convective 


regions at the core boundary) have been presented in Papics 
(2014 2015| ); Saio et al. (20151; Van Reeth & et. al. 


et al. 


(20151. None of them clearly exhibit a quadratic dependence 


on n, or a transition from a larger period spacing to a smaller 
one. This may provide evidence that semi-convective regions 
do not occupy large portions of low-mass stellar interiors, as 
found by ( Moore Garaud|2015 I. 

Although we have focused on g-modes, one may also 
wonder how a semi-convective staircase affects the propa¬ 
gation of p-modes. It is more difficult to obtain analytical 
results in this case, but it is still possible to apply the ma¬ 
trix transfer formalism used for g-modes. In so doing, we 
generically find that the leading order fractional correction 
to p-mode frequencies is 0[{d/HY\, which is much less than 
for g-modes when the wavelength is less than the scale height 
(i.e. A <C H). Moreover, unlike for g-modes, there is no lower 
cutoff frequency introduced by the finite size of the steps 
(though there is still a lower cutoff at N, the analog of the 
Brunt-Vaisala frequency introduced in equation §) . Both 
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of these reasons make g-modes much better candidates than 
p-modes for identifying regions of semi-convection in aster- 
oseismic power spectra. 

An effect which we have not included in our analysis, 
but which is likely to be important in astrophysical semi¬ 
convection is the forcing of modes by convective eddies. In 
particular, the presence of convection within a stair should 
preferentially excite modes with frequencies that are of order 
the eddy turnover frequency. Indeed, excitation of modes 
has been observed in simulations of layered semi-convection 
(e.g. [Wood et al.| ( |2013| )) and measuring the frequencies of 
such modes in simulations would be a useful check on the 
assumptions used to derive our analytical results. 
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APPENDIX A: JUMP CONDITIONS AT THE 
INTERFACES IN THE PRESENCE OF 
ROTATION 

We derive jump conditions at the interfaces between steps in 
the presence of a constant rotation rate. We start by writing 
down the linearized fluid equations in the Boussinesq ap¬ 
proximation in the corotating frame. The direction of strat¬ 
ification is along the 2 axis, and without loss of generality, 
we take fl — Q^x + fl^z. The primes denote differentiation 
with respect to z, and we have Fourier transformed in x, y, 
and t: 

0 = ikxSu + ikySv + 5w' (Al) 

iLoSu = - 2Q„5v 

P 

iujSv = _ 2flj:Sw + 2fl2(5it 

P 

■ r SP' or. r ^P 

lujow =-h zilxOV + g — 

P P 

. 'ySp dlnPp ' 

lUJ - = —ow -;- 

p dz 

Note that as in !|^ we have ignored the centrifugal force in 
the equations of motion, which assumes the rotation rate is 
far below breakup. 


Al Rotation Axis Along 2 

When the spin axis is along the direction of stratification 
{Qx = 0) we can set ky = 0 without loss of generality. The 
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set of equations ( |A1[ ) can be combined into a pair of first 
order partial differential equations, 


1 - 




Sw = 


SP = 


SP' 

iujp 

iujp 

1.2 


(A2) 


1 - ^ ) Sw', 


where N{z) is defined via equation 0 and the primes de¬ 
note differentiation with respect to The set of equations 
( |A2| | can be combined into a single second order partial dif¬ 
ferential equation for 5w. 

pk^ {ud — 5w = {ud — {p5w')'. (A3) 

Because equation( |A3[ ) is a second order differential 
equation, we need two jump conditions on Sw across the 
interface. Using the subscript “0” for the fluid above the 
interface and the subscript “1” for the fluid below the in¬ 
terface, the hrst of these conditions is Swo = Swi, which 
expresses the requirement that the upper and lower fluids 
stay in contact across the interface. 

To obtain the second jump condition, we integrate equa¬ 


tion (A31 across the interface to obtain 


-g—Sw = --(^^0 - <^^1) ■ 


(A4) 


Here, we have taken Ap = pi—po and assumed that Ap/p ^ 
1. The gA.plp term comes from integrating across the 
interface. 

The interface conditions (Isl are then modihed to read 


- B„+ie 

A„+ie“ -f = A„ + - B„), 

p 


(A5) 


where k is given by equation (56l. The dispersion relation 
(|57[) follows directly. 


smaller than g, under the Boussinesq assumption. Note that 
ujky/{k^ -|- fcy) < Loj ,,Jk'^ + ky <C Cs, where Cs is the sound 
speed. Since the scale height is H ~ c^/p, we can write: 


r\jx ' 


(A9) 


The last inequality follows from the fact that < fl^R 
g, which comes from our assumption that the centrifugal 
force is dynamically unimportant. Thus, equation ( |A8[ ) sim¬ 
plifies to 

,2 

(Swq — Sw'i). (AlO) 


OJ 


Ap, _ 

^ p kl + kl 


The interface conditions ([^ are then modified to read 
= A^-B,, (All) 


A D ^ — 

.^n-t-ie IJn + lG 


A I D ^ — _ A I D \ ^ ij^x T Ap ^ ^ T> \ 

An + ie -|- ±)„ + ie — An-\-rSn-\ - = — 7. - [An — Un), 

kuj^ P 


where k is given by equation ( 661 . The dispersion relation 


(671 follows directly. 


A2 Rotation Axis Perpendicular to z 

When the spin axis is perpendicular to the direction of strat¬ 


ification (Hj = 0), the set of equations (A11 can be combined 


into the following pair of first order partial differential equa¬ 
tions 


1 _ ^ _ 


4H 


^ SP' ^ 

Sw = - -h . ^ SP 

tup lUJ^p 


(A6) 


SP = - 


ip 


U2 I 

rvj. T rx,y 


(2flxkySw — LuSw') . 


The set of equations (A 61 can further be combined into a 


single second order partial differential equation for Sw: 

p {ud — N'^ — 4f2^) Sw = 

(—4fcyH^p — 2kyQ,xwp'^ Sw + LtP' {pSw')' 


1.2 I U2 


(A7) 


To obtain the jump condition for Sw', we integrate equa- 
intei 


tion (A71 across the interface to obtain 

,2 


2ky^\j,LjJ 
1.2 I U2 


+ 9 


p kx H” ^y 


(Swo — Sw'i) . (A8) 


We can simplify equation ( |A8[ ) further by showing that 
the hrst term in parentheses on the left hand side is much 
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